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Abstract 

We determine the location Ac of the mobility edge in the spectrum of the hermitian Wilson 
operator in pure-gauge ensembles with plaquette, Iwasaki, and DBW2 gauge actions. The results 
allow mapping a portion of the (quenched) Aoki phase diagram. We use Green function techniques 
to study the localized and extended modes. Where Ac > 0 we characterize the localized modes 
in terms of an average support length and an average localization length, the latter determined 
from the asymptotic decay rate of the mode density. We argue that, since the overlap operator is 
commonly constructed from the Wilson operator, its range is set by the value of A“^ for the Wilson 
operator. It follows from our numerical results that overlap simulations carried out with a cutoff 
of 1 GeV, even with improved gauge actions, could be afflicted by unphysical degrees of freedom 
as light as 250 MeV. 
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I. INTRODUCTION 


Domain-wall and overlap fermions reconcile chiral symmetry with the lattice, allow¬ 
ing for exact chiral symmetry at hnite lattice spacing in the euclidean path-integral 
formulation [1-5]. While chiral symmetry can be achieved for a range of non-zero bare 
coupling qq, problems arise if the bare coupling is too large. For domain-wall fermions, chi¬ 
ral symmetry cannot be maintained in the strong-coupling limit [6]. For overlap fermions, 
the built-in (modified) chiral symmetry is exact, but at strong coupling one loses either 
locality [7, 8] or control over the number of species [4, 9]. 

In any numerical simulation it is important to stay away from the dangerous regions of 
the phase diagram. The lattice Dirac operators of domain-wall fermions (DWF) and overlap 
fermions are both based on a Wilson operator with a negative, super-critical bare mass 
Locality and chirality in these formulations are controlled by the spectral properties of this 
underlying Wilson operator. 

The outstanding features of the super-critical Wilson operator are best illustrated in 
a theory with two dynamical flavors of Wilson fermions. Here the absence of a spectral 
gap in part of the phase diagram implies the existence of propagating, light degrees of 
freedom. Moreover, a non-zero spectral density for vanishing eigenvalue signals spontaneous 
symmetry breaking, as follows from the Banks-Casher relation [10]. For Wilson fermions 
there is no chiral symmetry to be broken, so the spontaneously broken symmetry is vectorial. 
As discovered by Aoki [11], the pions become massless if the bare Wilson-quark mass rriQ 
is lowered from positive values towards a critical value ruo = mdgo) < 0. For tuq < 
rndgo) the curvature of the effective potential for pions becomes negative at the origin; a 
pionic condensate forms which breaks spontaneously isospin and parity. This is the Aoki 
phase. Inside the Aoki phase the condensing pion is massive, while the other two pions are 
Goldstone bosons of the spontaneously broken isospin generators. See Fig. 1 for a schematic 
representation of the phase diagram. 

For DWF or overlap fermions one aims for a big gap, of order 1/a, in the spectrum of 
the Wilson operator. Simply speaking, a bigger gap in the spectrum of the Wilson operator 
improves both the chiral symmetry of DWF (at fixed, hnite extent of the hfth dimension) 
and the locality of overlap fermions. (As we will shortly see, in reality any such gap can be 
mostly, but not completely, devoid of eigenvalues.) Assuming that the phase diagram of the 
underlying Wilson operator remains qualitatively the same as in Fig. 1, this requires being 
inside one of the C phases. In practice, the rightmost C phase is used which, for go —0, 
coincides with the interval —2 < aniQ < 0. In this interval, one lattice domain-wall (or 
overlap) held gives rise to one quark in the continuum limit. 

When one studies the spectral properties of the Wilson operator as a kernel for DWF or 
overlap fermions, one is in fact considering a quenched Wilson-fermion theory, because the 
Boltzmann weight is derived from a diherent fermion operator. Usually, quenching means 
leaving out the fermion determinant altogether, but we will also consider the more general 
sense that the fermion determinant is that of DWF or the overlap operator. In this paper, we 
will study the spectrum of the Wilson operator for a variety of pure-gauge theories, but we 
will take our results as indicative of what might arise in a theory with dynamical fermions. 

The Goldstone theorem connects spontaneous symmetry breaking with the appearance of 


^ The super-critical region is — 8 < amo < 0. Outside of this region the Wilson operator cannot have zero 
eigenvalues. 
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FIG. 1: A schematic representation of the Aoki phase diagram. The lightly shaded area is the 
super-critical region. Phase B is the (massless) Aoki phase, while phases A and C are massive. 
Depending on the action in the theory, each of the thin “fingers” where the Aoki phase touches 
the line go = 0 could alternatively be replaced by a line of first-order phase transition [12]. Recent 
numerical results support the latter scenario in the case of dynamical Wilson fermions [13]. The 
phase diagram is symmetric under the replacement amo —> — (8 -|- amo), which can be undone by 
redefining the helds. Dotted lines indicate where we performed our numerical analysis. 


massless poles in correlation functions. Quenching the Wilson fermions, however, opens the 
door to a new dynamical possibility. Considering again a two-flavor theory, let us suppose 
that all the correlation functions of the two (quenched) Wilson flavors decay exponentially. 
(As we will see, this is indeed the case well inside the C phases.) In the quenched theory, 
this does not preclude a non-zero pion condensate: The isospin symmetry can be broken 
spontaneously without creating a Goldstone boson. It was shown in [7, 14] how this can be 
reconciled with the usual Ward-identity argument for the existence of a Goldstone pole. 

There is now solid numerical and semi-analytical evidence [15, 16] for the existence of zero 
modes of the Wilson operator throughout practically the entire super-critical region. The 
Banks-Casher relation again leads to a non-zero pion condensate. In particular, the pion 
condensate is non-zero for parameter values in the C phase that yield good-quality DWF 
and overlap-fermion simulations. Despite the condensate, all the Wilson-fermion correlation 
functions are short ranged. 

Ref. 7 provides a theoretical explanation of this situation. The low-lying eigenvectors 
of the hermitian Wilson operator may be either extended or (exponentially) localized. In 
the first case, the condensate must be accompanied by Goldstone bosons. In the second 
case, there is an alternative mechanism for saturating the relevant Ward identity, and all 
correlation functions can be (and, in fact, are) short-ranged. This gives the following physical 
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picture for the quenched Wilson-fermion phase diagram.^ In the entire super-critical region 
there is no gap in the spectrum of the Wilson operator, and the pion condensate is non-zero. 
For eigenvalues A above a certain mobility edge Ac, the eigenvectors of the Wilson operator 
are extended. If Ac > 0, eigenvalues |A| < Ac correspond to localized eigenmodes. In this 
case the pion condensate is non-zero, but there are still no long-range correlations. When 
Ac = 0, on the other hand, the condensate arises from extended eigenmodes, and there are 
Goldstone pions. 

The quenched Aoki phase is identihed with the region where Goldstone pions exist; that 
is, it is dehned by Ac = 0. With this dehnition, the quenched phase diagram could be 
qualitatively similar to that depicted in Fig. 1. Early numerical evidence supporting this 
quenched phase structure may be found in Ref. 18. The weak-coupling region may also be 
studied via an effective lagrangian [19]. 

As far as DWF and overlap fermions are concerned, the requirement of a gap in the 
Wilson spectrum should be replaced by the requirement that Ac > 0 [7]. In other words, 
one must work outside of the (quenched) Aoki phase, in one of the G phases. It is therefore 
important to map out the Aoki phase on any ensemble used for DWF or overlap-fermion 
numerical simulations. Furthermore, for practical reasons, one should not be too close to 
the Aoki phase. How close is “too close” depends on the underlying Boltzmann weight and 
on the construction of the fermion operator. We will discuss this very practical point at 
some length in our conclusions. 

In this paper we study the spectral properties of the Wilson operator via calculation of 
its resolvent and correlation functions derived from it. The theoretical framework 

developed in Ref. 7 is directly applicable, and guides us in the numerical implementation. 
We measure the spectral density as well as properties that characterize the shape and size of 
the localized eigenmodes. The resolvent gives us these quantities much more economically 
than would the direct study of the eigenvalues and eigenvectors of Hw- The correlation 
functions address the Ward identities and the Banks-Gasher relation directly. 

The resolvent allows simultaneous treatment of localized and extended modes. In any 
volume, the eigenvalues corresponding to localized modes are random. When the resolvent 
is averaged over the gauge ensemble, the single-conhguration spectral density, V~^ ~ 

\n), is smeared. Thus the ensemble-averaged spectral density is a continuous function of 
the eigenvalue for localized as well as extended modes. The essential physics of the localized 
modes lies not in their discreteness but in the compactness of their wave functions. 

Our measurements are carried out on pure-gauge ensembles (which is the usual meaning 
of quenching). We compare the spectral properties for three different pure-gauge actions: 
the standard plaquette action and the Iwasaki [20, 21] and DBW2 [22, 23] actions, two gauge 
actions motivated by renormalization-group considerations. These gauge actions have been 
used in quenched [24, 25] and dynamical [26] DWF simulations, and in quenched overlap 
simulations [27, 28]. 

This paper is organized as follows. In Sec. 2 we give basic dehnitions and derive the 
relevant Ward identities. In Sec. 3 we review the Banks-Gasher relation as well as the 
localization alternative to Goldstone’s theorem. A twisted-mass term [11, 29] provides the 
“magnetic held” that determines the direction of the pion condensate. Gareful study of the 
vanishing twisted-mass limit reveals that, if the low-lying eigenmodes are localized, the two- 
point function of the would-be Goldstone pions diverges linearly with the inverse twisted 


^ See also Ref. 17, in particular Fig. 2 therein. 
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mass. This enables the relevant Ward identity to be saturated without a Goldstone pole. 

We then turn to our numerical investigations, starting with the standard plaquette action 
for the gauge held. In Sec. 4 we present results for the simplest quantity, the spectral density. 
In Sec. 5 we dehne the localization length and use it to determine the mobility edge Ac for 
several points ( 5 ^ 0 , ^ 0 ) in the phase diagram. Extrapolations of Ac to zero allow us to map out 
a part of the boundary of the Aoki phase. We then proceed to detailed study of the localized 
modes. In Sec. 6 we extend the investigation to the Iwasaki and DBW2 gauge actions. We 
conclude in Sec. 7 with a discussion of the implications of our results for domain-wall and 
overlap fermions. The results of Sec. 5 yield several quantities that help locate regions of 
the phase diagram to be avoided in simulations. 

A concise account of this work, not including the improved gauge actions, has already 
been given [30]. 


II. DEFINITIONS 


A. Fermion action 


The Wilson-Dirac operator is dehned as 


D{mo) = 


I f iW + amo) 
a Gt 


(W + amo) ) ’ 


where 


9 ^ \^^x+fi,yUxfi 6x-(i^yUy^ 

^ y 


(Tn 


comes from the naive Dirac operator, and 

1 
2 


Wxy — 4:6xy ^x+jl,yUxy, + ^x-fi,yUyy 




( 2 . 1 ) 

( 2 . 2 ) 


(2.3) 


is the Wilson operator that breaks chiral symmetry while preventing species doubling. Here 
where at are the three Pauli matrices; we are using a chiral basis for the Dirac 
matrices, where 75 is diagonal. Uxy. is the SU(A) matrix representing the gauge held. We 
study the spectrum of the hermitian Wilson-Dirac operator. 


Hwijno) = D(mo) 75 - 

The corresponding eigenvalue equation (in a given gauge held) is 


(2.4) 


Hw{jno)^n — An'hn, 


(2.5) 


and we normalize the eigenvectors according to a^Yhx l'J'n(a^)P = 1 - 

Previous studies of the spectrum of Hw [15, 31, 32] were based on the calculation of 
individual eigenfunctions and eigenvalues. We hnd it more economical to calculate the 
Green function 

G{z) = [Hw{mo)-z]-\ ( 2 . 6 ) 

where z = X + imi, in order to extract information about the spectrum. G{X + imi) is well 
dehned in hnite volume provided mi 7 ^ 0. It has the spectral representation 


G{x,y,X + imi) = ^ 

n 


Xn — X — imi 


( 2 . 7 ) 
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B. Two flavors and twisted mass 


As mentioned in the Introduction, the spectral properties of Hw have profound effects 
on the realization of continuous symmetries when there is more than one flavor. Thus we 
will add an isospin index to the fermion held and consider the two-flavor theory dehned by 

Sp = — A)'^^ 

= - Ayg)^/’, (2.8) 

where tjj' = 75 '^. Spontaneous breaking of the flavor symmetry (and of parity) will be 
connected with the condensation of the “pion” held, 

P±{x) = i'lp{x)'-f^T±'4){x) 

7r3(a;) = i'4){x)'y^T‘i'4){x), (2.9) 

where t± = (ti -^ 1 x 2 ) The parameter A has been introduced into Eq. (2.8) in order to shift 
the focus from zero to nonzero eigenvalues of Hw- In order to control the isospin orientation 
of the condensate, we add to the action a “magnetic held” in the guise of a twisted-mass 
term, giving Anally 

Sp = il)[Hw - (A -f iT^mi)\il)' 

= - (A-f ir3mi)75]^/>. (2.10) 

rrii will be used as a regulator to avoid the singularities of G{z) along the real axis. 

C. SU(2) flavor symmetry and Ward identities 

For nil = 0, the fermion action (2.10) has a (vector) SU(2) flavor symmetry. For mi 7 ^ 0, 
the Ward identity of the broken symmetry is obtained by performing a local flavor transfor¬ 
mation, 

ipi^x) —>• ipi^x) + a{x)5+%l){x), ( 2 . 11 ) 

and similarly for ip{x), where 

5+'iIj{x) = ir+'^(a;), 5+ip{x) = —iip{x)T+. (2.12) 

We find for any operator O that 

7 + 2m, KU)O(s)) = '-H {s^o(y)). (2,13) 

(X 

Here the backward lattice derivative is defined by d*^f{x) = [f{x) — f{x — fi)]/a, and the 
vector current corresponding to Eq. (2.12) is 

\ \^{x)T+{'yf, - l)Uf,{x)ij{x + fl) + ^{x + + l)Ul{x)ij{x)\ . (2.14) 

While the notation (...) indicates an integration over both fermion and gauge fields, the 
Ward identity (2.13) in its various guises is in fact valid for each gauge configuration sepa¬ 
rately. 
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We define the “pion” two-point function 


V{x,y) = {Ti+{x)Ti-{y)) (2.15) 

as well as 

(a:, y) = (j+(a;)7r_ (?/)). (2.16) 

[Note that this is the usual pion only when A = mi = 0; see Eq. (2.10).] If we take 
Oiij) = Ti-iy) in Eq. (2.13), we find the Ward identity 

'9;)r^(a^, y) + 2 mir(a;, I/) = ^ (7^3(1/)) • (2.17) 

Applying a Eourier transform, viz. 

f(p) = (2-18) 

h(p) = (2.19) 

we derive the momentum-space Ward identity, 

- ^(1 - e-“7’-)f ^(p) + 2mif (p) = (tts) . (2 .20) 

III. GOLDSTONE’S THEOREM AND LOCALIZATION 

The Ward identity (2.20) is valid for arbitrary A and mi. In a quenched theory, however, 
despite the Goldstone Theorem, (tts) 7 ^ 0 does not necessarily imply the existence of a 
massless pole in r^(p) in the limit mi —0. Let us recall [7, 14] how this comes about. 


A. Condensate and Banks Casher relation 

The volume-averaged pion condensate 

(tts) = {a^/V) Y, ( 773 ( 2 :)), (3.1) 

X 

in the two-flavor theory (2.10) can be expressed in terms of the Green function G(A ± imi). 
We will denote the expectation value in a given gauge held by (.. .)u- Then 

(i*'Y ~ —Tr [G(A-|-imi) — G(A — fmi)] 

where we have used the spectral representation (2.7). Averaging this over the gauge held 
gives the translation-invariant result 

(3.3) 
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(3.4) 


where p(A) is the eigenvalue density defined by 

^ n 


In the limit mi 0, we obtain 

(tts) = 27rp(A). (3.5) 

This is a generalized Banks-Casher relation; the original relation [10] is Eq. (3.5) at A = 0. 


B. Localization as an alternative to Goldstone’s theorem 

Naively taking the limit mi —0 in Eq. (2.20) gives 

(>^ 3 ) = 

® fl 

- (3.6) 

where the second line is the approximate form for ap -C 1. As we shall see shortly, Eq. (3.6) 
is sometimes false, but it contains the Goldstone Theorem; r^(p) must have a massless pole 
for any A such that (tts) ^ 0. By the generalized Banks-Casher relation (3.5), this happens 
whenever p(A) 7 ^ 0. Apart from T^{x,y), one expects long-range power-law decay also for 
other correlation functions, including in particular r{x,y). 

In the physics of disordered systems [33] it is well known that the eigenmodes of a hamil- 
tonian in a random background divide into two classes: extended and localized. In fact, the 
spectrum splits into bands, each band containing only eigenmodes of one type.^ A point in 
the spectrum separating an extended band from a localized band is a mobility edge. r^(a;, y) 
will exhibit a power-law decay when the A lies in an extended band, while for a localized 
band it will decay exponentially. We expect that the same basic separation applies to Hw 
as well. 

If p{\) comes from localized eigenmodes and r^(p) has no pole at zero momentum, what 
has become of the Goldstone Theorem? In other words, how is the Ward identity (2.20) 
satisfied? The way out of this conundrum is the following. In the limit mi 0, r(p) 
diverges as 1/mi for a range of values of p that includes the point p = 0. The limiting value 
of mir{p) is finite. For p —0, we arrive at an alternative to Eq. (3.6), 

lim 2mif(0) = lim (tts) . (3.7) 

mi^O mi^O 


C. Divergence of the pion two-point function 

Let us consider further the 1/mi divergence in r(p). The (finite-volume) spectral repre¬ 
sentation of the charged-pion two-point function is 




An- — A -|- imi 


<-(l/)^n+(l/) 


A„+ — A — imi 


(3,8) 


^ There seems to be no rigorous proof of this fact except in one dimension [33]. 



where terms with the subscript are associated with the propagator for the corresponding 
quark flavor. As explained in Sec. 3 of Ref. 7, a 1/mi divergence may arise only from the 
terms with = n~, so that 

In analogy with Eq. (3.4), we define the eigenmode-density correlation function, 

n,ix,y) = |d/„(a;)nvl/„(l/)|25(A - A^), (3.10) 

n 

and its Fonrier transform. 


*a(p) = Ry - A„)), (3.11) 

where 

^n(p) = a^El^n(a;)|V^’A (3.12) 

X 

As mi —A 0, these may be interpreted as the contribntion to T{x,y), or to its Fonrier 
transform, of eigenmodes with eigenvalue A. Repeating the analysis leading from Eq. (3.2) 
to Eq. (3.5) we hnd, for mi 1, 

r(i,j/) = ^2AM)+0(1). (3.13) 

mi 

Observe that lZ\{x,y) is the ensemble average (3.10) of a quantity that is strictly positive. 
Also, since the eigenmodes are normalized, we have that 

’^a(0) = ^1:'Ra(x,p)=p(A), (3.14) 

and so p(A) ^ 0 means that 7l\{x,y) must be nonzero for at least one pair of valnes {x,y). 
Consequently, in any finite volnme, p(A) 7 ^ 0 implies the existence of a 1/mi divergence in 
the coordinate-space two-point fnnction. [This must be true for at least one valne of {x,y), 
bnt is expected to hold for practically every {x,y).] This resnlt is valid for the generalized 
quenched theory dehned by Eq. (2.10) for any Boltzmann weight; the only assnmption we 
used is that the Boltzmann weight does not depend on A and mi. 

[An unquenched theory with fermion action (2.10) wonld include dei{Hw — A — iT^mi) in 
the Boltzmann weight. This determinant will snppress eigenvalues A„ ~ A, and the spectral 
density measured by Eq. (3.5) will be zero in any hnite volume in the limit mi —a 0.] 

We next consider the inhnite-volnme limit. The asymptotic behavior of an exponentially 
localized eigenmode is 

|T^(a;)|^ - exp ^ ^nl j ^ (315) 

where is the localization length and x^ is the center of the localized eigenmode. Eq. (3.15) 
is valid only at distances that are large compared to the size of the region containing most 
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of the eigenmode’s density. In principle, nothing forbids the occurrence of eigenmodes with 
a very short localization length < a. 

An extended eigenmode is one that does not satisfy Eq. (3.15) for any hnite Evidently, 
truly extended eigenmodes exist only in inhnite volume. In hnite volume, the clear-cut 
identihcation of an eigenmode as localized demands that is exponentially small on 

most of the lattice. In later sections we will give a more quantitative criterion. 

In inhnite volume, we expect the Fourier transform of the eigenmode density (3.15) to 
have the following small-p behavior; 

„ipxl 

^ 11 2 / 2 - ( 3 . 16 ) 

1 + pHi 

The region 1 will rehect only the exponentially decaying envelope (3.15) of the 

eigenmode density but not the short-distance huctuations. The overall normalization is set 
by Hn{0) = 1 for a normalized eigenmode. Substituting the ansatz (3.16) into Eq. (3.11) 
and going to small p we hnd, in analogy with Eq. (3.13), 


f(p) 


’rp(O) 

mi 


l + 0(p"7")] +0(1), 


(3.17) 


where I is some average localization length for the eigenmodes with eigenvalue A. [This 
extends Eq. (3.7) to nonzero p.] From Eq. (2.20) we conclude that 

f^(p) ~ ip^p{0)f. (3.18) 


Thus there is no Goldstone pole when the eigenmodes with eigenvalue A are exponentially 
localized.^ 

When the eigenmodes at the given A are extended, the transition from Eq. (3.8) to 
Eq. (3.9) is not justihed if the limit mi —0 is taken after the inhnite-volume limit, because 
of interference effects between eigenmodes with inhnitesimally close eigenvalues. In this case 
we expect Eq. (3.6) to be valid, indicating a Goldstone pole. 


IV. WILSON PLAQUETTE ACTION: SPECTRAL DENSITY 

We now turn to our numerical investigation. We begin with quenched ensembles gener¬ 
ated with the Wilson plaquette action for the SU(3) gauge theory, 

A=|^ReTr(l-t/,^,). (4.1) 

O X 

^J.<V 

We began calculating at /? = 6.0, which is usually taken to correspond to the lattice scale 
a~^ ~ 2 GeV, and at^ mo = —1.5, between the Aoki “hngers” that point to the P = oo 
{do = 0) axis at mo = 0 and —2. Then we moved downward in P towards the Aoki phase, 
calculating at /? = 5.85, 5.7 (where a~^ ~ 1 GeV), 5.6, 5.5, and 5.4. We will show that the 
Aoki phase is entered just below P = 5.6 (see Fig. 1). 


^ See also Sec. 4 of Ref. 7. 

^ For the remainder of the paper we rescale amo —>■ mo, giving the bare Wilson mass in lattice units. 
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Returning to (3 = 5.7, we moved “sideways” by changing mo to —2.0 and —2.4. The 
latter turns out to be very close to or in the second Aoki finger. We put these choices of 
(/5, mo) into the context of other work in Sec. VII. 

For each value of f3 we generated 120 uncorrelated gauge configurations (except where 
otherwise noted) on a lattice of 16"^ sites, using the MILC pure-gauge overrelaxation code. 
For a given (3, results for all values of mo. A, and mi were calculated on the same ensemble; 
thus correlations had to be taken into account in all fits and statistical analysis. 

We shall illustrate our methods by discussing in detail the analysis for (3 = 5.85 and 
mo = —1.5. Resnlts for other valnes of {(3, mo) will be summarized in the tables. 


A. Spectral density from Banks—Casher relation 

From the Banks-Casher relation (3.5) and Eq. (3.2) we have 

2np{X) = lim (tts) = ^ lim (ImTr G(A + imi)). (4.2) 

mi^O y mi^O 

Thus the volume- and ensemble-averaged Green function, extrapolated to mi = 0, gives 
p{X) directly. Of course mi must be kept nonzero for actnal calculation in order for G to be 
bounded. 

The spectral snm (3.2) shows how to do the extrapolation. For any gange configuration 
U we consider the snm 

y__ (4 3) 

y(A.-A)2 + mf 

The snmmand tends to a (5-fnnction as mi —0, but before the limit is taken it has a finite 
width eqnal to mi. The given configuration will make a contribntion of 0(l/mi) if Hw 
has an eigenmode whose eigenvalne A^ satisfies |A„ — A| ^ mi; these contribntions, snmmed 
over configurations U, will add np to a finite limit as mi —> 0. On the other hand, all the 
eigenmodes that are far from A, with a|A„ — A| =0(1), will make a contribntion of 0(mi). 
This indicates a linear extrapolation, 

= Co + cimi, (4.4) 

Zn 

where co,i will depend on A. Then cq is an estimate for p(A). 

We calculated Tr G using a single random sonrce per gange configuration. Averages were 
obtained for np to seven values of my. 0.01, 0.02, ..., 0.07. The upper graph in Fig. 2 shows 
the linear extrapolation for two values of A. For A = 0.3 the fit works well. For A = 0., 
on the other hand, the extrapolated intercept is very small and the precision attained is 
inadeqnate. 

The problem is that, when p(A) is too small, there are few eigenmodes (for the volume 
we use) within the broadened h-function in Eq. (4.3). Then the snm is dominated by the 
0(mi) contributions of the more distant eigenmodes, with large finctuations. Making mi 
even smaller will suppress these contributions, but it will be self-defeating because even 
fewer modes will lie within the broadened h-fnnction, increasing the finctnations in their 
contribntion as well. 
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FIG. 2: Top graph: the pion condensate (27r)“^(7r3) as a function of mi (open symbols), and its 
linear extrapolation to yield the spectral density p(A) at mi = 0, for A = 0.3 (upper curve) and 
A = 0 (lower curve). Bottom graph: cubic extrapolation of improved estimator for A = 0 (circles), 
compared to the result of the linear fit above (square). All data are for Wilson action at /3 = 5.85; 
mo = —1.5. 

B. Improved estimator 

A solution lies in using an improved estimator for the spectral density, one that suppresses 
the contribution of distant eigenmodes and thus approaches the mi —0 limit faster than 
linearly. We introduce the (dimensionless) differential operator J dehned by 

Jf = (4.5) 

orrii \mi J 
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for any /(mi). It is designed to remove the linear term in /(mi). Applying it to the pion 
condensate, we have 

2 

J (tts) = — (^Im Tr G — mi Re Tr , (4.6) 

with the spectral representation 

The mi —0 limit again yields 27rp(A) since / dxm\/{x^ + m\Y = 7r/2. 

The operator J indeed removes the leading, linear contribution of the eigenmodes with 
a|An — A| = 0(1). We can see from Eq. (4.7) that the contribution of these eigenmodes is 
now proportional to m\. In view of this, we attempt a cubic extrapolation, 

= Co + csm?. (4.8) 

111 

The lower graph in Fig. 2 shows the new extrapolation, again at A = 0. Linear and quadratic 
terms are in fact unnecessary to attain an excellent £t to the data points; this is the final jus- 
tihcation of the model (4.8). Most important, the precision in the extrapolation is improved 
by a factor of ten. 

Using the improved estimator (4.6) requires a second inversion for each random source. 
The number of CG iterations needed for the second inversion can reach twice the number 
for the first inversion. The cost of the improved estimator is thus up to (roughly) three 
times the original cost. Were we to invest the additional computer time instead in increased 
statistics using the simple estimator (4.2), the anticipated reduction in the statistical error 
would not come anywhere close to the factor of ten achieved with the improved estimator. 


( 4 . 7 ) 


C. Spectral density from the two-point fnnction 

Since the two-point function T{x,y) is the eigenmode-density correlator [see Eq. (3.9)], 
it is a rich source of information about spectral properties. The first application we discuss 
is an alternative method to calculate the spectral density p(A). 

We project r(a;, y) to zero spatial momentum and calculate the time-correlation function, 

rw = Trllr(0,x;t,y), (4.9) 

^3 x,y 

where V 3 is the three-volume. This calculation requires a random source on time slice 0, as 
well as a random sink for each time slice t] again we use one set of random sources per gauge 
configuration. If we sum r{t) over all t, we can use translation invariance and Eq. (3.9) to 
show that 

lim ami ^ r(t) = 7 rp(A). (4-10) 

This corresponds to setting p = 0 in the finite-volume Ward identity (2.20). [Note the 
similarity to Eq. (3.14).] We verified Eq. (4.10) for two values of the bare coupling {j3 = 6.0 
and 5.7) by calculating both the condensate [Eq. (4.2)] and the two-point function. In all 
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subsequent measurements we did not calculate the condensate separately since it gives only 
the spectral density, while r(t) yields this and much more. 

The extrapolation of mir{t) to mi = 0 suffers from the same fluctuation problems as that 
of (tts) when the spectral density is small. The solution again lies in an improved estimator, 
but now for the eigenmode-density correlation function 7l{x, y, A) [see Eq. (3.10)]. We dehne 
the “improved” two-point function, 

rimp(a:, y) = Tr {g{x, y)Qiy, x )), (4.11) 


where 


Q{x,y) = a‘^ y~] G(x, z] \ + imi)G{z, y, \ — imi) 

Z 

^ ^ T^(a;)T),(|/) 

„ (A„ - A)2 -hmf ’ 

and project it to zero spatial momentum, 

a® 

rimp(t) = — ^rimp(0,x;f,y). 

V?, 


3 x,y 


(4.12) 


(4.13) 


Instead of Eq. (3.9) one has 


whence 

TT 

\im mlTimp{x,y) = -nx{x,y). 

mi^O 2 

Equation (3.14) then yields the spectral density [cf. Eq. (4.10)], 

'jj- 

lim am? ^ri,^p(t) = -p(A). 


(4.14) 

(4.15) 


(4.16) 


Results obtained using T(t) were extrapolated to mi = 0 linearly in mi, while those obtained 
using riinp(t) were extrapolated® linearly in mf. 

We also tried to improve the estimator for the eigenmode density by calculating j7T(a;, y) 
[cf. Eq. (4.5)]. This has the advantage that the application of JT” to all terms in the Ward 
identity (2.20) generates a new identity. As it turns out, while some reduction in the 
statistical error was achieved this way, the improvement was not signihcant. The reason is 
that the application of J' to T{x,y) would be effective in suppressing the contribution of 
eigenmodes with a|A„± — A| ~ 1 only if the approximation Eq. (3.9) is valid. As it turns 
out, for the mi values we used, Eq. (3.9) is not a good approximation to Eq. (3.8) when 
the spectral density is small. The correlator rimp(a;, y) does a better job in suppressing the 
contribution of eigenmodes with a|A„± — A| ~ 1. 


The form of Eq. (4.14) suggests an term as well, but we found it to be unnecessary for a good fit. 
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D. Numerical results 


We present our results for the spectral density for the Wilson plaquette action in Tables 
1-8 and in Figs. 3 and 4. 

At each value of f3 and mo, we performed our calculations for a range of values of A 
starting from zero, increasing in steps of 0.1. In the ranges studied (the maximal A was 
0.4 to 0.6) the spectral density typically increases by one to three orders of magnitude. For 
the ensemble generated by the plaquette action, p(A) shows no remarkable behavior as A 
passes Ac, the mobility edge (to be determined below). In all cases, however, the mobility 
edge is encountered for p(Ac) ~ 0.1. There is a steady rise in the spectral density as we 
decrease f3 (Fig. 3), which bespeaks an increased disorder in the gauge held. 

A hnal note on improvement: For the larger values of A it is possible to achieve acceptable 
precision without using the improved estimator. Since the number of CG iterations required 
grows rapidly with the spectral density, we limited the use of improved estimators to those 
cases where results obtained without improvement were poor. 


TABLE I: Spectral properties for (3 = 6.0, mo = —1.5. The mobility edge is at Ac 0.41, marked 
by the horizontal line in the table. The quantities V, Ig, and R, which characterize the localized 
modes, are defined in Sec. V. 


A 

P 

h 

V 

Is 

R 

0.0“* 

0.00011(2) 

0.61(1) 

0.000038(8) 

2.9(8) 

31. 

0.1“^ 

0.00024(2) 

0.61(1) 

0.00010(1) 

2.4(3) 

25. 

0.2“ 

0.0016(1) 

0.70(1) 

0.00053(6) 

3.0(4) 

16. 

0.3“ 

0.0214(9) 

1.16(4) 

0.0058(4) 

3.7(3) 

8. 

0.4“ 

0.128(3) 

3.8(6) 

0.025(2) 

5.1(4) 

5.3 

0.5 

0.267(7) 

oo 

0.033(5) 

8.(1) 

4.4 

0.6 

0.390(9) 

oo 

0.042(9) 

9.(2) 

4.0 

“Results from improved estimator. 





^Used 1200 gauge configurations for measurements with rrii 

= 0.01, 0.02, 0.03. 




TABLE 11: Spectral properties for (3 = 5.85, mo = 

— 1.5. The mobility edge is at Ac — ( 

3.32. 

A 

P 

k 

V 

Is 

R 


0.0011(1) 

0.64(1) 

0.00036(5) 

3.1(5) 

17. 

0.1“ 

0.0019(1) 

0.71(1) 

0.00054(7) 

3.5(5) 

15. 

0.2“ 

0.0088(4) 

0.92(2) 

0.0024(2) 

3.7(3) 

10. 

0.3“ 

0.056(1) 

2.2(2) 

0.0132(8) 

4.2(3) 

6.5 

0.4 

0.168(7) 

oo 

0.036(6) 

4.7(8) 

4.9 

0.5 

0.27(1) 

oo 

0.023(7) 

12.(4) 

4.4 

0.6 

0.39(2) 

oo 

0.055(9) 

7.(1) 

4.0 


“Results from improved estimator. 
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TABLE III: Spectral properties for (3 = 5.7, mo 


1.5. The mobility edge is at Ac — 0.25. 


A 

P 

h 

V 

Is 

R 

0.0“ 

0.0089(3) 

1.04(3) 

0.0023(2) 

3.9(4) 

10. 

0.1“ 

0.0126(7) 

1.29(5) 

0.0029(3) 

4.3(5) 

9.4 

0.2“ 

0.039(1) 

2.3(1) 

0.0097(5) 

4.0(2) 

7.1 

0.3“ 

0.106(2) 

oo 

0.026(1) 

4.1(2) 

5.5 

0.4“ 

0.199(6) 

oo 

0.036(4) 

5.5(6) 

4.7 

0.5“ 

0.308(8) 

oo 

0.052(4) 

5.9(5) 

4.2 


“Results from improved estimator. 


TABLE IV: Spectral properties for (5 = 5.6, mo = —1.5. The mobility edge is at Ac = 0.14(2). 
(We give the statistical error in Ac here because it is larger than in other cases. All suffer from 
systematic error in the extrapolation.) 


A 

P 

k 

V 

Is 

R 


0.027(2) 

2.1(1) 

- 

- 

7.8 

0.1 

0.034(3) 

4.(1) 

- 

- 

7.4 

0.2 

0.072(4) 

oo 

0.012(3) 

6.(2) 

6.1 

0.3 

0.135(7) 

oo 

0.016(5) 

8.(3) 

5.2 

0.4 

0.219(8) 

CO 

0.017(6) 

13.(5) 

4.6 


TABLE V: Spectral properties for f3 = 5.5, mo = —1.5. The mobility edge is at zero. 


A 

P 

k 

V 

Is 

R 

o 

o 

0.057(1) 

oo 

0.0112(6) 

5.1(3) 

6.5 

0.1“ 

0.073(1) 

oo 

0.013(1) 

5.6(4) 

6.1 

0.2 

0.116(5) 

oo 

0.015(4) 

8.(2) 

5.4 

0.3 

0.172(7) 

oo 

0.019(4) 

9.(2) 

4.9 

0.4 

0.25(1) 

oo 

0.029(6) 

9.(2) 

4.5 


“Results from improved estimator. 


TABLE VI: Spectral properties for f3 = 5.4, mo = —1.5. The mobility edge is at zero. 


A 

P 

k 

V 

Is 

R 

0.0 

0.066(3) 

oo 

- 

- 

6.2 

0.1 

0.089(5) 

oo 

0.008(3) 

11.(4) 

5.8 

0.2 

0.148(6) 

oo 

0.022(4) 

7.(1) 

5.1 

0.3 

0.194(7) 

oo 

0.018(4) 

11.(3) 

4.8 

0.4 

0.27(1) 

oo 

0.028(6) 

10.(2) 

4.4 
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TABLE VII: Spectral properties for /3 = 5.7, mo = 

—2.0. The mobility edge is at Ac ^ 

0.21. 

A 

p 

k 

V 

Is 

R 


0.014(2) 

1.14(4) 

0.005(2) 

3.(1) 

9.1 


0.022(2) 

1.37(5) 

0.005(2) 

4.(2) 

8.2 


0.13(2) 

3.8(5) 

0.03(1) 

4.(2) 

5.3 


0.24(1) 

oo 

0.028(4) 

9.(1) 

4.5 


0.41(2) 

oo 

0.05(1) 

8.(2) 

4.0 


TABLE VIII: Spectral properties for /3 = 5.7, mo = 

—2.4. The mobility edge is very close to zero. 

A 

P 

k 

V 

Is 

R 


0.043(2) 

4.(1) 

0.006(2) 

7.(2) 

6.9 

0.1 

0.112(5) 

oo 

0.014(3) 

8.(2) 

5.5 

0.2 

0.234(8) 

oo 

0.026(6) 

9.(2) 

4.5 

0.3 

0.40(1) 

oo 

0.052(6) 

8.(1) 

4.0 

0.4 

0.54(1) 

oo 

0.05(1) 

11.(2) 

3.7 


V. WILSON PLAQUETTE ACTION: LOCALIZATION PROPERTIES 

A central goal of the work presented here is the determination of the mobility edge Ac at 
various places in the (/5, mo) phase diagram. Equally interesting is the shape and density of 
the localized eigenmodes below Ac (where Ac > 0). 

The two-point function r{x,y) contains complete information about the eigenmodes of 
Hy/-, and it is our primary tool both in determining Ac and in studying the localized modes. 
When dealing with localized eigenmodes, we aim for parameter values where the single 
spectral sum (3.9) provides a good approximation to the exact expression (3.8). When mi is 
small enough, r(a;, y) then reduces to the eigenmode-density correlator IZx {x,y) [Eq. (3.10)]. 
Our analysis is based on this feature. The approximation (3.9) ceases to hold when the 
eigenmodes become too dense or too extended, that is, when they interfere in the double 
sum. This occurs when A is above or too close to the mobility edge. We will develop a 
criterion to establish the consistency of our analysis of the localized eigenmodes. 

We begin by giving a precise dehnition of the localization length /^(A). Its divergence 
marks the mobility edge Ac(/5,mo). By dehnition, Ac = 0 marks the Aoki phase. We then 
turn to study other properties of the localized eigenmodes outside the Aoki phase, where 
Ac > 0. 

A. The localization length 

We have loosely dehned the localization length of an individual eigenmode in Eq. (3.15). 
We can use /„, nebulous as it is, to motivate the dehnition of an average localization length 
that can be obtained from from the large-f behavior of r(t). Once we reach this dehnition, 
a precise dehnition of In will be superhuous. 
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FIG. 3: Eigenvalue density p(A) for (bottom to top) f5 = 6.0, 5.85, 5.7, 5.6, 5.5, 5.4. Mobility edges 
Ac(/3) (when Ac > 0) are indicated by vertical bars. Wilson action; niQ = —1.5. 


We begin by introducing the restricted spectral density pz(A) that contains the contribu¬ 
tions of eigenmodes with localization length < I only. It is given by 

= W-U). (5.1) 

and its derivative, the differential spectral density, is given by 

^ = 1? ( E5(A„ - min - I)). (5.2) 

We also introduce a probability distribution V\{1) for the localization length I of eigenmodes 
with given eigenvalue A by writing p((A) = p{\)V\{l). This distribution is normalized 
because, in hnite volume, 

[ dl p'i{X) = p{X) and hence [ dlV\{l) = 1, (5.3) 

Jo Jo 

where L is the (largest) linear size of the lattice. Below, this upper limit will be implicit. 

[In inhnite volume there can be truly extended eigenmodes. Since we expect the eigen¬ 
modes at a given A to be either all extended or all localized, we have correspondingly either 
V\{1) = 0 for all hnite I, or dlV\{l) = 1.] 

The decay rate of individual localized eigenmodes can be related to the decay rate of r(t) 
for mi ^ 0 [7]. Using Eq. (3.15) in Eq. (3.9) gives 
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FIG. 4: Eigenvalue density p(A) for /3 = 5.7 and (bottom to top) mo = —1.5, -2.0, -2.4. Mobility 
edges are indicated by vertical bars. For rriQ = —2.4 the mobility edge is probably very close to 
A = 0. Wilson action. 


miT{x,y) ~ E exp 


X - x'^\ + \y - x[ 


nil 




Averaging over the position of gives 


miT{x,y) ~ (^exp| 


k - y\' 

Irt 


(A„ - A)2 -h ml 


mi 


- A)2 ml 


2 / ’ 


(5.4) 


(5.5) 


since the average is dominated by locations near the straight line connecting x and y. 
Hence 

\x-y\\ 


lim miVix^y) ~ 7rp(A) / dlV\(l) exp 
mi^O J ' 


I 


where we have used Eq. (5.2). Equation (4.9) then gives 

lim miT(t) ~ 7ip{X) [ dlV\{l) expf— 
mi^O J \ 


(5.6) 


(5.7) 


Thus r(t) decays exponentially, as a weighted average of exponentials with all localization 
lengths. Motivated by this result and ignoring power corrections, we hnally dehne the 
average localization length as 


hW = 


p{xy 


(5.8) 
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FIG. 5: Squared decay rate of the two-point function r(t) vs. mi (open symbols) and its linear 
extrapolation to mi = 0. Top to bottom: A = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6. The mobility edge is at 
A ~ 0.32. All data are for Wilson action at /3 = 5.85; mo = —1.5. 


where /r(A) is the extrapolation to mi = 0 of the decay rate M of the two-point function, 

r(f) ~ exp[—M(A, mi)f], t a. (5.9) 

Equation (5.7) actually suggests a slightly different dehnition of /r(A), based on the ex¬ 
trapolation of the correlation function as a whole, 

lim mir(t) ~ exp[—p(A)d, t S> a. (5.10) 

mi^O 

In practice, the extraction of /i(A) using Eq. (5.10) turns out to be extremely noisy, but 
compatible with the results obtained using Eq. (5.9).^ We therefore used the dehnition 
based on Eq. (5.9) throughout our calculations. 

We extrapolate to mi = 0 by htting 

M^(A, mi) =/i^(A)-f Q;(A)mi. (5-11) 


^ Some further details are provided in Ref. 30. 
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See below for an explanation of this choice. We illustrate this extrapolation in Fig. 5. The 
tables list the results for li = for comparison with other characteristic lengths.® 

The derivation leading to Eq. (5.7) relies on Eq. (3.9), where interference effects present 
in Eq. (3.8) have been dropped. This is certainly true when the exponentially localized 
eigenmodes are isolated. A precise dehnition of what this means will be given later. For 
now it suffices to say that, intuitively, the eigenmodes with eigenvalue near A are isolated if 
every such eigenmode is supported in a different part of the lattice. 

As A increases, the spectral density grows and so does the localization length of individual 
eigenmodes. Eventually, interference sets in, meaning that the approximation of Eq. (3.9) 
is not valid for the mi values we use. Interference will cause F(t) to decay faster than 
the decay rate of individual eigenmodes. The dehnition Ig = will then systematically 
underestimate the true (average) localization length, as extracted perhaps by calculating 
individual eigenmodes and matching them to Eq. (3.15). If we are to use Eq. (5.9) to 
calculate the localization length, we must verify that the eigenmodes are indeed isolated. 
We return to this issue below. 


B. The mobility edge 

As we increase A, we reach a critical value Ac where /u(Ac) = 0 (see Fig. 5). What is 
the physical signihcance of Ac? According to the preceding discussion, /i“^(A) provides a 
reasonable average value for the localization lengths of individual eigenmodes if there is 
no interference, while it underestimates the average when interference sets in. Either way, 
/i(A) = 0 implies that the average localization length of individual eigenmodes is inhnite; 
the eigenmodes have become extended. The point Ac therefore marks the mobility edge. 

Would a different procedure yield a different value for the mobility edge? One might 
consider, for instance, an alternative determination based on the numerical calculation of 
individual eigenmodes. We believe that such ambiguity can only be an artifact of hnite- 
size effects. Such effects are perforce signihcant when /i(A) < 1/L, since the notion of 
extended eigenmodes is truly meaningful only in inhnite volume. Thus some disagreement 
between different determinations of the mobility edge is to be expected in any hnite volume. 
Nonetheless, we expect that all methods will converge to the same value in the inhnite- 
volume limit [33]. We are only interested here in obtaining an overall picture, which does 
not justify the resources that would be needed for, say, a hnite-size scaling analysis. 

We estimate Ac simply by a linear extrapolation of fj,{\) to zero, using the last two non¬ 
zero values we have found. The justihcation for this procedure is that, for A near Ac, we 
expect that /i(A) will exhibit the characteristic behavior of a (continuous) phase transition. 
See Figs. 6 and 7. 

Above the mobility edge we expect the 1/mi divergence of T{x,y) to disappear in the 
inhnite-volume limit; hence Eq. (3.6) is valid, and F^(p) should have a Goldstone pole. We 
expect that the physics of these Goldstone poles to be governed by some effective chiral 
lagrangian. Borrowing from the familiar physics of ordinary Goldstone bosons, we expect 
oc mi in the presence of the symmetry-breaking held mi. This is the motivation for the 
term linear in mi in the ht function (5.11). Below the mobility edge does not vanish for 


® We calculated the localization length from r(t), without improvement. Where the tables indicate im¬ 
provement, it applies to p and V only. 
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FIG. 6: The inverse squared localization length jj? vs. A. Top to bottom: (3 = 6.0, 5.85, 5.7, 5.6. 
Intercepts at the A axis are extrapolated from the last two points in each set, and provide estimates 
for the mobility edge Ac. Data are for Wilson action; thq = —1.5. 


mi —> 0, and therefore it is immaterial whether we extrapolate M or to mi = 0. Put 
together, these considerations suggest that Eq. (5.11) is an appropriate form both below 
and above the mobility edge. 

Above the mobility edge, the linear extrapolation often gives a small negative result 
for /i^. We have neglected to include possible logarithms which, in the context of a chiral 
lagrangian, would occur at next-to-leading order in mi. These logarithms, together with 
finite-size corrections, should move the extrapolated value to zero. 

Our linear extrapolation to hnd the zero of fi{X) = suffers from an uncertainty that 

stems from the interference to which we alluded above. On further analysis (see Sec. VE 
below) we will conclude that, for all (/?, mo) where Ac > 0, our estimate of is unreliable for 
the largest value of A that is below the mobility edge Ac- (This uncertainly does not occur 
for lower values of A.) In fact, the calculated at this value of A is likely to be too high; 
its true value could even be zero. In other words, this value of A could lie above the mobility 
edge. This means that the last measured point on each curve in Figs. 6 and 7 is unreliable. 
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FIG. 7: Same as Fig. 6, at /3 = 5.7. Top to bottom: mo = —1.5, —2.0, —2.4. For mo = —2.4 the 
mobility edge is probably very close to A = 0. Wilson action. 


C. Entering the Aoki phase 

Our estimates for the mobility edge Ac(/3, mo) are collected in Table IX. Consider first 
the mo = —1.5 results. For reference, we include the free-theory limit (/3 = cx)) [7], where Ac 
coincides with the gap of the free Hw- At /3 = 6.0, Ac is still close to its free-held value. As 
we move towards strong coupling, the curve Ac(/3) steepens before reaching zero somewhere 
between /3 = 5.6 and j3 = 5.5. As discussed in Sec. Ill, where Ac = 0 the two-point function 
r(t) has a zero-momentum pole even at A = 0. The Goldstone Theorem is valid, the Ward 
identity (3.6) is satished, and the (quenched) Wilson-fermion theory possesses a Goldstone 
boson associated with the spontaneous breakdown of its SU(2) havor symmetry. This is the 
Aoki phase, and we put its boundary between /3 = 5.6 and (3 = 5.5. 

Table IX also shows results at the other two mo values for (3 = 5.7. We hnd only a small 
change between mg = —1.5 and mo = —2.0, while the mo = —2.4 result suggests that here, 
too, one is near or within the Aoki phase. 

The implications of these results for DWF and overlap fermions are discussed in Sec. VII. 
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TABLE IX: Mobility edge Ac(/3, mo). Where no error is shown, the statistical error is less than one 
in the last digit. 


p 

mo 

Ac 

oo 

-1.5 

1/2 

6.0 

-1.5 

0.41 

5.85 

-1.5 

0.32 

5.7 

-1.5 

0.25 

5.6 

-1.5 

0.14(2) 

5.5 

-1.5 

0.0 

5.4 

-1.5 

0.0 

5.7 

-2.0 

0.21 

5.7 

-2.4 

ss 0 


D. Participation number and the support length 


Now we turn to a more detailed description of the localized eigenmodes, for which |A| < Ac- 
We begin by dehning a measure of the size of the support of a localized eigenmode. In the 
sequel we use this measure to sharpen the notion of isolated localized eigenmodes, and verify 
the consistency of our approach. 

We dehne [31, 33] the participation number P of a normalized eigenmode 'h(a;) in d 
dimensions by 

p-i=a‘'i:(|>l-(nPL (5.12) 

X 

The physical meaning is easy to see. Suppose that the magnitude of the eigenmode density 
is |\l/(a;)p ~ l//f over a region whose linear size is Is- Then P^P ^ is a measnre of the 
linear size of the snpport of \l/(a;). We similarly dehne a generalized participation nnmber 
by 

X ( X l>l'(x„)|f. (5,13) 

[Eq. (5.12) corresponds to d' = 0.] In the case considered we wonld have 


p—^ l^d—d^^l^d^—d'^2 _ l^' 


(5.14) 


In the special case d' = d — 1, we have Pd,d-i ~ h- 

We are thus motivated to dehne ls{n) = to be the linear size of the snpport of 

the eigenmode. Snbstitnting in Eq. (3.9) and nsing translation invariance together with 
Eqs. (4.9) and (5.14) we have 


Dehning 


we have 


mir(t = 0) 


e(? ls{n) (A„-A)2 + m?) + ^^"''^- 

(5.15) 

P(A) = 7i~^ lim mir{t = 0), 

mi—^O 

(5.16) 

®(A) = p(A) JdlViiDj, 

(5.17) 
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where Vl{l) is the probability distribution for Isiji). On the basis of Eq. (5.17) we dehne 
the average support length as 

aA) = ^. (5.18) 

While we expect the average localization length and the average support length to be 
quantities of similar magnitude, there is no reason why they should be the same. 

We may also obtain 'D{\) using the improved estimator (4.13). Again, results obtained 
using T(t) were extrapolated to zero linearly in mi, while those obtained using riinp(t) were 
extrapolated linearly in mf. 

We list our results for Ig in the tables. Observe that below the mobility edge turns out 
to be always larger than The physical signihcance of these results is discussed below. 


E. Separation distance of localized eigenmodes 

We now return to the notion of isolated localized eigenmodes mentioned above. The use 
of Eq. (3.9) depends on neglecting interference among the modes appearing in Eq. (3.8), 
which is only justihed when mi is sufficiently small. We must check whether our smallest 
values of mi are indeed sufficiently small. 

Spectral sums as in Eq. (3.9) show that mi is the resolution with which we detect eigen¬ 
modes near A. Let us compare the average support and localization lengths Ig and Ig to the 
mean distance between eigenmodes detected at this resolution. Our smallest value of mi is 
0.01. The number of eigenmodes with eigenvalues near A that we detect for a typical gauge 
conhguration is thus N{\) ~ O.OlEp(A). Hence we dehne 

R{X) = [V/N{X)]^/^ ~ [0.01p(A)]“^/^ (5.19) 

as a measure of the average separation between eigenmodes detected at this resolution. 
Results for R{X) are shown in the tables. (We do not quote errors but they can be easily 
worked out from the p(A) data.) If Ig R the eigenmodes are isolated, and correlation 
functions rehect properties of individual localized eigenmodes, with no interference.® Thus, 
our method for the extraction of the average localization and support lengths is valid. 

We can also directly estimate the overlap among any two eigenmodes to see whether they 
interfere. This overlap satishes the inequality 

I (^i|^2) P < l^i(^)ll^2(a:)|)'. (5.20) 

X 

We assume Ai ~ A 2 ~ A. The “edge” of the support of 4/2 closest to the center of Ti is at 
an average distance of around R — lg/2. Using Eq. (3.15) this gives 

I (T 1 IT 2 ) P < exp[-(R - lg/2)/k]. (5.21) 

A review of the tables shows the following.^® Leaving out the last value of A just below 
the mobility edge we hnd that for all other A < Ac one has Ig < R/2 and Ig < R/5. Hence 


® Here it is significant that Ig turns out to be smaller than Ig. 

This applies also to results obtained with improved gauge actions, to be discussed in Sec. VI below. We 
ignore those cases where data for L are not available. 
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FIG. 8; Coefficient of the 1/mi-divergence in for (3 = 5.7, mo = —1.5, A = 0.0. 

I (\l/i|\l/ 2 ) P ~ exp[—(3i?/4)/(i?/5)] = exp(—15/4) ~ 0.024. This demonstrates that, except 
for the last value of A just below Ac, all the quantities we measured indeed reflect properties 
of individual localized eigenmodes, with no interference. 

Above Ac, as well as (in most cases) for the A value just below Ac, we have R ^ Ig ^ 
This means that interference is not a negligible effect, and the values of Ig and Ig no longer 
reflect the properties of individual eigenmodes. 

F. The 1/mi divergence of the two-point function: p/ 0 

Finally, we return to the Ward identity Eq. (2.20). As can be seen from Eq. (3.17), the 
localization alternative to Goldstone’s theorem requires that the 1/mi divergence persist for 
a range of momenta, and that its coefficient depend smoothly on p. 

In order to confirm this, we calculated the Fourier transform of F(f), 

f(cUn) = Xlcos(cUnt)F(t) (5.22) 

t 

where oon = 27rn/16, and extrapolated mir(ci;„) linearly to rrii = 0. [The justification for 
the linear extrapolation is the same as for Eq. (4.4). We included data for mi = 0.005 in the 
fit.] The calculation, which provides a consistency check, was done for f3 = 5.7, mo = —1.5, 
A = 0, where the higher spectral density (compared io j3 = 6.0) makes it easier to obtain 
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good statistics; the mobility edge here is still far from zero. We plot the results in Fig. 8. 
The cu-dependence of the l/rrii divergence is indeed smooth, as indicated in Eq. (3.17) for 
small uj. 

For comparison, we repeated the calculation for A = 0.5, which lies above Ac. The 
extrapolation of = 0) to mi = 0 is again straightforward, as it must be since this 

gives p(A) according to Eq. (4.10). Doing the same with r(ci;„ 7 ^ 0), however, leads to a 
huge showing that Eq. (3.17) is inapplicable. This highlights the qualitative difference 
between |A| < Ac and |A| > Ac. 


VI. IMPROVED GAUGE ACTIONS 

As an alternative to the Wilson plaquette action, we also studied pure-gauge ensembles 
generated by the Iwasaki [20] and DBW2 [22] actions, which have been used in DWF [24-26] 
and overlap [27, 28] simulations. We set the Wilson mass to mo = —1.8 (see Sec. VII for 
explanation of this choice). For each action, we studied two values of the lattice spacing, 
a~^ ~ 2 GeV and a~^ ~ 1 GeV. The corresponding bare couplings are (3 = 2.6 and 2.2 
respectively for the Iwasaki action [21]. For the DBW2 action they are (3 = 1.04 [25] and 0.79 
respectively; the latter was determined by interpolation using /3-values from Ref. 23. 


TABLE X: Spectral properties for Iwasaki gauge action, (3 = 2.6, mo = —1.8. The mobility edge 
is at Ac 0.38. Number of configurations ranges from 600 at A = 0. to 100 at A = 0.5. 

A 

P 

k 

V 

Is 

R 

o 

o 

7.(3) X 10-'^ 

0.72(1) 

— 

— 

no. 

0.1“ 

1.8(5) X 10-5 

0.73(1) 

4.(1) X 10-6 

4.(2) 

50. 

0.2“ 

2.3(9) X 10-4 

0.80(1) 

1.3(7) X 10-4 

2.(1) 

25. 

0.3“ 

9.(2) X 10-4 

1.20(2) 

2.7(8) X 10-4 

3.(1) 

18. 

0.4“ 

0.111(2) 

oo 

- 

- 

5.5 

0.5“ 

0.277(5) 

oo 

0.040(4) 

7.(7) 

4.3 

“Results from improved estimator. 


TABLE XL Spectral properties for Iwasaki gauge action, f3 = 2.2, mo = — 
is at Ac cx 0.27. Measurements on 200 gauge configurations. 

1.8. The mobility edge 

A 

P 

k 

V 

Is 

R 

o 

o 

0.0063(2) 

0.80(1) 

0.0021(1) 

3.0(2) 

11. 

0.1“ 

0.0095(3) 

0.91(2) 

0.0032(2) 

3.0(2) 

10. 

0.2 

0.034(2) 

1.40(5) 

0.010(2) 

3.4(7) 

7.4 

0.3 

0.120(5) 

oo 

0.025(4) 

5.(1) 

5.4 

0.4 

0.27(1) 

oo 

0.017(7) 

16.(10) 

4.4 


“Results from improved estimator. 
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We present our results in Tables 10-13 and Figs. 9 and 10. A comparison of the (inverse 
squared) localization length for the three gauge actions at 2 GeV is shown in Fig. 9. The 
Iwasaki and DBW2 localization lengths are almost the same, an indication of the good 
scaling properties of the two actions. The mobility edge is almost the same for all three 
actions. Our results at 1 GeV are compared in Fig. 10. There is a gradual decrease in 
the localization length from Wilson to Iwasaki to DBW2. Differences in the value of the 
mobility edge, while bigger than in the 2 GeV case, remain small. 

As can be seen from Tables 10-13, the main difference among the three actions is a 
dramatic reduction in the low-A spectral density as we go from the Wilson action to the 
Iwasaki and then to the DBW2 gauge actions. What distinguishes the three actions is the 
coefficient of the 1x2 rectangle, which is zero for the standard plaquette action and in a 
common parametrization [24, 25] is c = —0.331 for the Iwasaki action and c = —1.4069 
for the DBW2 action. We see that the decrease in the low-A spectral density is correlated 
with a more negative rectangle coefficient c. The rectangle term in the action suppresses 
the small-size dislocations [22, 25] responsible for the existence of the near-zero modes. 

At (or just above) the mobility edge the spectral density remains of the same order as for 
the plaquette action: p(Ac) ~ 0.1. This means that, unlike the plaquette action, the Iwasaki 
and DBW2 spectral densities rise steeply as A approaches the mobility edge from below. At 
2 GeV, a comparison of the spectral-density values just above and below the mobility edge 
shows a jump by two orders of magnitude for both improved actions. 


TABLE XII: Spectral properties for DBW2 gauge action, (3 = 1.04, niQ = —1.8. The mobility edge 
is at Ac ^ 0.39. Between 200 and 400 gauge configurations were used. 

A 

P 

h Is 

R 

0.0“ 

- 

0.74(1) 

- 

0.1“ 

2.(1) X 10“^ 

0.75(1) 

150. 

0.2“ 

6.(2) X 10-6 

0.82(1) 

65. 

0.3“^ 

3.4(7) X 10-^ 

1.2(2) 

23. 

0.4 

0.058(2) 

oo — — 

6. 

“Results from improved estimator. 

^Extrapolation with Eq. (4.8) failed. We added terms linear and quadratic in mi. 



TABLE XIIL Spectral properties for DBW2 gauge action, (3 = 0.79, uiq = 
edge is at Ac — 0.32. Measurements on 200 gauge configurations. 

— 1.8. The mobility 

A 

P 

k 

V 

Is 

R 

0.0“ 

0.0013(1) 

0.70(1) 

0.00041(5) 

3.2(5) 

17. 

0.1“ 

0.0023(2) 

0.73(1) 

0.0009(1) 

2.6(4) 

14. 

0 2“fe 

0.0064(4) 

0.91(1) 

0.0024(2) 

2.7(3) 

11. 

0.3“ 

0.058(1) 

2.1(1) 

0.012(1) 

4.8(4) 

6.4 

0.4 

0.22(1) 

oo 

0.030(5) 

7.(1) 

4.6 


“Results from improved estimator. 

^Extrapolation with Eq. (4.8) failed. We added an term. 
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FIG. 9: The inverse squared localization length vs. A, as in Fig. 6. Comparison of three gauge 
actions at a~^ ~ 2 GeV: Wilson action at /3 = 6.0, with mo = —1.5 (circles); Iwasaki action at 
(3 = 2.6, with mo = —1.8 (squares); DBW2 action at (3 = 1.04, with mo = —1.8 (diamonds). 

VII. DISCUSSION 

Our numerical results bear, first, on the phase structure of the super-critical Wilson 
operator in various ensembles; and, second, on the overlap and DWF theories based on the 
Wilson operator. We examine each of these in turn. 


A. The Aoki phase diagram 

The Aoki phase diagram has dictated our list of (/3, mo) values for calculation, and our 
results in turn yield information about the diagram. We can now discuss in more detail the 
choice of (/3, mo) values we have made. 

There is considerable evidence (see e.g. Ref. 15) that, for plaquette actions with (3 = 5.7 
to 6.0, the low-A spectral density of the Wilson operator varies slowly over the range —2.0 < 
mo < —1.5. In view of this, we have chosen to carry out most of our calculations with the 
plaquette action for the single value mo = —1.5. We have tested this insensitivity to mo by 
considering as well (3 = 5.7, mo = —2.0. We see in Fig. 7 that the localization lengths li for 
all A are almost unchanged as mo changes from —1.5 to —2.0; hence, the mobility edge is 
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FIG. 10: As in Fig. 9, but at a~^ ~ 1 GeV: Wilson action at /3 = 5.7, with mo = —1.5 (circles); 
Iwasaki action at /3 = 2.2, with mo = —1.8 (squares); DBW2 action at /3 = 0.79, with mo = —1.8 
(diamonds). 


unchanged as well. The rough equality extends to the support length and to the spectral 
density for A sufficiently below the mobility edge, which means in practice A < 0.1 (see 
Tables III and VII as well as Fig. 4). 

Old results [18] on the Aoki phase diagram indicate that the point (/d = 6.0, mo = —2.4) 
lies in the second “finger” of the Aoki phase. Our results show that the point (5.7, —2.4) is 
indeed at a boundary of this finger or inside it. 

Looking ahead to DWF, our results indicate that, for f3 > 6.0, one can use any value 
of mo in the range [—2.0,—1.5]. (The same applies to the overlap construction.) Mean- 
field arguments suggest that the optimal mo (related to the “domain-wall height” M via 
M = —mo) for DWF simulations is roughly mo ~ —1 -|- rndgo) [34]. The commonly used 
value at a~^ ~ 2 GeV is M = —mo = 1.8. This is consistent with the mean-field estimate 
and lies comfortably in the middle of the range of mo we considered. In view of our evidence 
for the insensitivity of the key spectral properties over this range, we have chosen the value 
mo = —1.8 also for the improved-action calculations. 

Our scan of (3 values for mo = —1.5 places the boundary of the Aoki phase of the plaquette 
action between (3 = 5.5 and [3 = 5.6 (cf. [18]). 
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B. Implications for overlap fermions 


The overlap operator is constructed explicitly in terms of the Wilson operator Hw, 

Oov = 1 - 

Our knowledge of the spectrum of is thus of immediate import to overlap calculations; 
in particular, it is critical to understanding the locality of the operator. The long-distance 
tail of -Dov has two components [7, 8], and both must be kept under control. One component 
comes from the near-zero modes; its range is the localization length /^(A), where A will be 
dehned shortly. The other component comes from modes near the mobility edge, and its 
range is the inverse of Ac. 

We denote by Ov the range of the overlap operator: For \x — y\ 3> a one has \Do^{x,y)\ ~ 
exp(—lx —|/|//ov)- In Ref. 7 we discussed the different roles of the extended and the localized 
modes in determining /qv, and arrived at an estimate for the contribution of the localized 
spectrum. Here we revisit the issue in the light of our numerical hndings. 

In the treatment of the localized eigenmodes the essential assumption was the absence of 
interference between different modes. This should be true, in particular, for the low-lying 
modes that make up the mutually isolated subspace. These are all the localized eigenmodes 
in the interval [—A, A] where 0 < A < Ac is chosen such that for any two eigenvalues Am, An G 
[—A, A], the corresponding localized eigenmodes are isolated (on average) in the sense of 
Sec. VE. The contribution of the mutually isolated subspace to the ensemble-averaged 
overlap operator is then estimated to be^^ 


(|L)ov(a:,|/)|)|;,|<3; ^ y^_dAp(A)exp • (7-2) 

Since both p(A) and /^(A) are monotonically increasing, the integral is dominated by its 
upper limit, viz. 

(|T)ov(a:,|/)|)|;,|<); ^ Ap(A)exp • (7-3) 

This is the hrst of the two pieces determining /qv 

We may arrive at a numerical estimate as follows. Consider for dehniteness our results 
for the Wilson gauge action at (3 = 6.0, shown in Table 1. A conservative guess for the edge 
of the mutually isolated spectrum is A ~ 0.2. Since p(0.2) is of the order of 10“^, we hnd 
that the integrated spectral density J^idXp{X) is of the order of 10“^. The total number 
of eigenvalues in the interval [—A, A] is given by the volume V times this integral, and so 
the mean space-time separation R between any two eigenmodes in this spectral intervaE^ is 
R ~ 10. Using similar considerations to Sec. VE, this corresponds to | P = 0(0.01), 


This estimate complements the discussion of Ref. 8. The analytical techniques of Ref. 8 are inapplicable 
to the localized spectrum because it becomes dense in the infinite-volume limit. 

The introduction of A improves on the discussion of Ref. 7, where the rather arbitrary cutoff Ac/2 was 
imposed on the integral [see in particular Eq. (6.3) therein]. 

Note that R measures the average separation between all modes in the interval [—A, A], whereas R{X) 
defined in Sec. VE measures the average separation between modes in an 0 (toi) neighborhood around A. 
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which is certainly a small overlap. All the modes in this spectral interval are, therefore, 
mutually isolated. 

Evaluating Eq. (7.3) for A = 0.2 gives 


(l^ov(a:,|/)|)|;,|<3, ~ 10 ^exp 

This semi-quantitative estimate highlights two points concerning A. First, the small prefactor 
signifies a large value of R. The latter must be appreciably larger than the support length 
Is in order for the modes to be mutually isolated, and A must be chosen to make it so. 
Moreover, pushing A to a somewhat higher value would also increase le{X) and thus modify 
the coefficient in the exponent. Equation (7.4) still provides a rough idea of the contribution 
of the mutually isolated localized modes. 

The modes with |An| > A are of two types: the localized eigenmodes in the interval 
A < I An I < Ac, and the extended modes above Ac. The rise in the spectral density near and 
above Ac leads us to expect that the resulting contribution to /qv will be controlled by Ac 
itself. For these modes the analysis of Ref. 8 applies, and a comparison of our results to 
theirs supports this conjecture. For f3 = 6.0 and mo = —1 — s = —1.4 and —1.6, Table 1 
in Ref. 8 reports ~ 0.49 and 0.45, respectively. This can be compared to our result for 
the mobility edge at mo = —1.5, namely, Ac — 0.41. Given the differences in the details, the 
close agreement may be a numerical coincidence. Still, this suggests that the contribution 
of the rest of the spectrum may indeed be dominated by the extended modes just above the 
mobility edge, and that 


(|T>ov(a:,i/)|)|;,|>X Cexp(-Ac|a; -yl), (7.5) 

where C = 0(1). This is the second contribution to the tail of Oov 

As explained above, the prefactor in Eq. (7.4) is bound to be small. Comparing Eqs. (7.4) 
and (7.5) suggests that /qv is dominated by the mobility edge and not by the near-zero modes. 
This is consistent with the results of Ref. 35. There it is found that changing the gauge 
action from Wilson to Iwasaki or to DBW2 at quenched ~ 2 GeV has little effect on 
lav- Indeed, while the spectral density of the localized modes is a sensitive function of the 
gauge action, the mobility edge itself varies little among the three gauge actions (see Tables 
10-13 and Fig. 9). Of course, if > AjT^ the small contribution of the near-zero modes 
will dominate the asymptotic tail of Hov at large distances. For all cases we study here, 
however, the opposite inequality holds. 

The overlap kernel’s decay rate can be interpreted as the mass of unphysical degrees 
of freedom. In practice, these degrees of freedom can be uncomfortably light. Assuming 
/oy Ri our data allow us to estimate their mass for the pure-gauge ensembles. When the 
cutoff is a~^ ~ 2 GeV, we find l~-^ ~ 0.4 x 2 GeV = 800 MeV for all three gauge actions. 
This is presumably a high enough scale to qualify as part of the discretization errors. For 
a~^ ~ 1 GeV, on the other hand, we obtain ~ 0.25 x 1 GeV = 250 MeV for the Wilson 
action, 270 MeV for Iwasaki, and 320 MeV for DBW2. This is an alarmingly low scale for 
unphysical degrees of freedom.^® 


This rough equality among the actions is in agreement with the above mentioned results of Ref. 35. 
Note in particular the 1 GeV overlap simulations of Ref. 27. 
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C. Implications for domain-wall fermions 


Our discussion of DWF will be less detailed for two reasons. First, the theoretical back¬ 
ground has already been discussed in detail [7, 36-38]. Also, the relevant “hamiltonian” for 
DWF is not H\y itself but the logarithm of the hfth-dimension transfer matrix, a different 
(though closely related) operator. 

DWF achieve exact chirality when the number of sites L 5 in the hfth dimension tends 
to inhnity. DWF simulations are performed at finite L 5 , typically in the range of 10-20. 
The main question is what is the size of chiral symmetry violations due to the finiteness 
of L 5 . A quantitative measure of these violations is provided by the residual mass rures 
[24, 25, 38], which is the small additive correction to the quark mass determined from the 
PCAC relation [ 2 ]. (Alternatively, m^es can be determined from the extrapolation of the pion 
mass to the chiral limit.) Here, too, there are two terms that can be ascribed to extended 
and localized modes. 


C-2 

nires ~ Cl exp(-AcT5) + 

-^5 


(7.6) 


where Ci = 0 ( 1 ) comes from the extended modes near the mobility edge and 


dXp{X)^p{l/L,) (7.7) 

J- 1 /L 5 

comes from the low-lying, localized modes. Because of the rapidly growing spectral density, 
the localized modes’ contribution is dominated by modes with |A,i| ~ I/T 5 ; we ignore a 
power-law correction to the extended modes’ contribution.^® The tildes indicate spectral 
quantities of the new “hamiltonian,” 

^ = -log(T2)/(2a5), (7.8) 


where T(a 5 ) is the transfer matrix for hopping in the hfth direction and 05 is the corre¬ 
sponding lattice spacing (conventional DWF have 05 = 1). 

The “hamiltonians” and H share identical zero modes [2]. Thus we expect that for 
the near-zero modes, the spectral density (and other properties) of H is fairly close to that 
of Hw- Replacing p(A) by p(A) in Eq. (7.7) should yield a reasonably good approximation. 
On the other hand, further away from A = 0 the spectra of H and Hw need not be equal. 
In particular the mobility edges could be quite different. Nonetheless, again because of the 
identity of the zero modes, the mobility edges of H and Hw reach zero simultaneously. Thus 
the same Aoki phase defines the forbidden region for both overlap and DWF. 

As noted, the numerical results of this paper are relevant for determining the near-zero 
modes’ contribution to rUres, but separate calculations would be needed to determine Ac, 
the mobility edge of H, which governs the extended modes’ contribution to rures- One can 


We thank N. Christ for a discussion of Eq. (7.7). This formula, derived from Appendix C.2 of Ref. 36, 
provides a slightly better estimate of the contribution of the localized modes than that given in Ref. 37. 
Equation (7.6) corrects the discussion of the near-zero modes contribution to Wres given in Ref. 7. In 
particular, Eq. (6.12) therein as well as the first term on the right-hand side of Eq. (7.1) there are erroneous. 
The limits of integration in Eq. (7.7) here come of demanding that exp(—AT 5 ) = 0(1) for modes included 
in the integral. 
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alternatively determine these quantities by calculating the residual mass as a function of L 5 
[24, 25, 38] and htting to Eq. (7.6). 

Equation (7.6) shows that, unlike overlap fermions, the physics of hnite-Ls DWF simu¬ 
lation and, in particular, the value of m^es, are sensitive functions of the near-zero modes of 
H (or of H\\r). In quenched DWF simulations, using the Iwasaki or DBW2 actions allows 
reaching negligibly small values of m^es [24, 25]. For overlap simulations, on the other hand, 
these gauge actions are advantageous for a technical reason: overlap simulations require 
an exact treatment (within numerical precision) of all the Wilson eigenvalues in a certain 
interval [—5,(5], and reducing the number of modes in this interval speeds up the simula¬ 
tion. We note that there exist versions of DWF, in particular the so-called Mobius fermions 
[39], where the near-zero modes’ contribution to mres decreases much faster with L 5 (see 
also [40-42]). 


D. Summary 

The continuum limit of lattice QCD with either DWF or overlap fermions can be taken 
while letting mo —1. Since all correlation lengths associated with the Wilson operator 
itself remain hnite in lattice units, there is little doubt that the continuum limit is correct, 
and that no unphysical excitations can survive it. Issues addressed in this paper have to do 
with MC simulations at hnite lattice spacing. 

Using Green function techniques, we have determined the mobility edge of the hermitian 
Wilson operator for a number of pure-gauge ensembles with plaquette, Iwasaki and DBW2 
gauge actions. Our results allow mapping a portion of the (quenched) Aoki phase diagram. 
Where the mobility edge is non-zero, we have also characterized the localized spectrum in 
terms of an average support length and an average localization length, the latter determined 
from the asymptotic decay rate of the mode density. 

Our results are of direct relevance to the overlap operator. For the near-zero modes, or, 
more precisely, for the mutually-isolated subspace of the localized modes, we found that the 
localization length is consistently smaller than the support length. We have also found that 
twice the localization length is not bigger than the inverse mobility edge. Together, these 
hndings imply that the near-zero modes play little role in setting the range of the overlap 
operator. 

We argue that the range of the overlap operator is set by, and is roughly equal to, the 
inverse mobility edge. Our results for the mobility edge suggest that it is fairly safe to 
perform (quenched) overlap simulations at a~^ = 2 GeV. The same is not true when the 
cutoff is 1 GeV; for all three gauge actions we hnd that the standard overlap operator is 
likely to contain unphysical quark-like degrees of freedom as light as 250-300 MeV. This casts 
serious doubt on the validity of 1 GeV overlap simulations. In general, the determination of 
the range of the overlap operator is an important test that must be carried out for any new 
overlap simulation. 

Glosely related to the mobility edge is the mass of the lowest pseudoscalar excitation 
of the Wilson fermion action—Wilson’s original would-be pion. It, too, is a non-physical 
excitation where the overlap operator is concerned. Where we demand that couplings be 
chosen such that the mobility edge is far from zero, the same can be said of the mass of the 
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Wilson pion.^^ 

Obtaining the corresponding information for DWF will require the study of a different 
“hamiltonian,” the logarithm of the hfth-dimension transfer matrix, H. In particular, it 
is important to study the range of the effective four-dimensional operator D^s obtained 
by integrating out the hve-dimensional bulk modes and the pseudofermions [40], and to 
determine how the range of Deff is affected by if as a function of L 5 and 05 . The near-zero 
spectrum of H is similar to that of the Wilson operator, and we are thus able to conhrm 
the picture that the near-zero modes (of either Hw or H) make a major contribution to the 
residual mass. 

Last, we note that numerical simulations with dynamical DWF [26] (or dynamical overlap 
fermions) require the largest available computers. Not only production runs, but also ex¬ 
ploratory runs are very expensive. The process of closing in on an optimal set of simulation 
parameters can greatly beneht from the (low-cost) determination of the quantities studied 
here—the mobility edge and the characteristics of the localized modes. 
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